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Abstract 



We propose a nuclear mechanism which can induce an atomic electric 
dipole moment (EDM). Parity and time reversal violating (P, T-odd) nuclear 
forces generate a mixed P, T-odd nuclear polarizability (3ik (defined by an 
energy shift U = —PikEiH^, E is electric field and H magnetic field). The 
interaction of atomic electrons with produces an atomic EDM. We per- 
formed an analytical calculation of the P, T-odd nuclear polarizability and es- 
timated the value for the induced atomic EDM. The measurements of atomic 
EDMs can provide information about P, T-odd nuclear forces and test models 
of CP-violation. 
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I. CALCULATION OF P, T-ODD NUCLEAR POL ARIZ ABILITY 



Time (T) reversal violation in the interactions of elementary particles has only been 
indirectly observed in the CP-violating decay of the neutral K-mesons. The observation 
of a P, T-odd multipole (for example, the dominant electric dipole moment (EDM)) of, for 
instance, an elementary particle or atom would be direct evidence for T-odd interactions. 
There are several schemes for incorporating CP- violation into unifying theories. These 
schemes make very different predictions for the value of EDMs of particles and atoms. 
While there has never been an unambiguous detection of these EDMs, upper limits on their 
values provide stringent tests of competing models of CP- violation. 

In this paper we propose a mechanism which can induce an atomic EDM: the interaction 
of a P, T-odd perturbed nucleus with the electric and magnetic fields of an external electron. 

The interaction of external electric E and magnetic H fields with operators of nuclear 
electric and magnetic dipole moments, d and p, gives rise to an energy shift, U = —fiikEiHk, 
where fiik is the mixed P, T-odd nuclear polarizability . Parity and time reversal violating 
(P, T-odd) nuclear forces generate this polarizability, 

p ^-2^2 (jgl^j^jjrj^^ (l) 

n E n — E 

where ipo and i[) n are P, T-odd perturbed ground and excited nuclear states, respectively. We 
consider a heavy spherical nucleus with a single unpaired nucleon. The P, T-odd interaction 
between the non-relativistic nucleon and the nuclear core is given by (see, e.g. [|J) 

G is the Fermi weak interaction constant, rj is a constant which characterizes the strength 
of the P, T-odd nucleon-nucleon interaction (limits on 77 have been obtained from atomic 
and molecular ||] experiments), a is the nuclear spin operator, p(r) is the density of 
the nuclear core, and m is the proton mass. The P, T-odd perturbed nuclear wavefunctions 
resulting from the interaction Hpt, in an approximation where the nuclear density p and 
potential U coincide, po(r) = U(r)p (0)/U(0), have the form [|J 

4=(l-^-V)Vn, (3) 

where £ = —rj 2 ^ m and ip n is an unperturbed wavefunction. 

In lowest order in the P, T-odd interaction the polarizability (Eq. [[]) is 

Pik = ~n< ErWol^ • V,d;]|V>n)(^n|/iA#o) + (^o|^|^n) (^n\ [<? ' V,/ifc]|^o)} (4) 

n l^n — 

= /4 1} + 4 2) - (5) 

In the non-relativistic approximation, the only contribution to the first term of this equation 
arises due to transitions between fine-structure doublets, for example py 2 and p 3 / 2 , d 3 / 2 and 
d-5/2 (because the only non-zero matrix elements of p are between these states). In the 
second term, states ip n must have opposite parity to tp (because the matrix element of 
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d = er cannot mix states of the same parity). The largest contributions to this term arise 
due to transitions to the next shell. Therefore the energy denominator of the second term 
is much larger than that of the first. 

The nuclear electric dipole operator is defined as d = e^r, where is the effective 
charge for the nucleon N = n,p, e p = (N/A)e, e n = —(Z/A)e, which appears due to the 
recoil effect. The nuclear magnetic dipole operator is defined as /x = fj,^r(gi\ + g s s), where 
fiN is the nuclear magneton, g s and gi are the spin and orbital g-factors; for the proton 
g s = 5.586, gi — 1 while for the neutron g s = —3.826, g\ = 0. 

The first term of Eq. [| can be expressed as 

Pik = e N /J-N^(9s - gi) — (^oWi\iJo>)(ipo>Wk\ipo), (6) 
where Uf s = E > — E is the energy difference between the fine structure components ipo and 
Let us now consider the second term of Eq. [3], f3\ k . This can be written as 

Pik = ^2~B — ^r{9s(ipo\d i \ipn}('ipn\Wj^k}^j\ipo} +2g l (ip \d i \^ n )(ilj n \a j [V j ,l k \\ilJ }}. (7) 

Using the commutation relations [o^Cfc] = 2iejk m u m and \Vj,h] = ifjfcmVm, we obtain 

Pik - Zl ^^NKg s - gijejmk 2^ B B • ^> 

n &n — 

This term can be further simplified by replacing the operator V by the commutation relation 
V = m [r, H], where H = p 2 /2m + V(r) is the single-particle Hamiltonian. Neglecting the 
spin-orbit interaction, that is setting [H,cr] = 0, and using closure, J2 n iVViXVVil = 1> Eq- H 
becomes 

pf k ] = -2i£e N /i N m(g s - 0i)ejy m (^oK^m|V'o)- (9) 

So, to first order in the P, T-odd interaction, we can write the nuclear polarizability as 

Pik = £e N fx N (g s - gi){(l/ujf s )(4)o\ai\ip ,)(tJj 0/ \a k \tJj ) - i2m(0o|?"i(<7 x r) k \ip )}. (10) 

The P, T-odd nuclear polarizability can be expressed in terms of scalar and tensor com- 
ponents, 

Pik = kkPs + [hh + hh - ^1(1 + 1)]/%, (11) 

I is the nuclear total angular momentum. We are interested in finding the scalar and tensor 
P, T-odd polarizabilities, P s and p t . Let us start with the scalar contribution. Now, from 
Eq. [11], we see that P s = Pu/3. We can see from Eq. [K^ that there is no scalar contribution 
from the second term, since rj(o~ x r)j = 0. We can therefore write the scalar term as 

Ps = £e N fi N (g s - gi)(^—)(ipoWi\ipo')(ipo'Wi\ipo)- (12) 

6UJ is 
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Evaluating this expression using spherical spinors, with the ground state wavefunction ipo 
corresponding to / = / + 1/2 and its fine-structure partner ip / corresponding to I' = I — 1/2, 
I is the orbital angular momentum, we obtain 

1 21 — 1 

P* = &nhn($8 - gi)( — )( ot )■ ( 13 ) 



Let us now consider the tensor polarizability f3 t . Setting i, k = z in Eq. [10] and taking 
the maximum projection, I z = I = I + 1/2, the first term is zero, since with maximum 
projection I' z — I' — I — 1/2 gives zero contribution. The second term is also zero in this 
case, since a x and a y cannot have diagonal matrix elements. Using Eq. [11] with j3 zz = we 
obtain 



Inserting Eqs. O and [T3 into Eq. O, we obtain for the nuclear P, T-odd polarizability 



= &NV N (g s - 9i)(—)(j) v\ ^ S ik - hlih + hh - h ik l(l + l)]}. (15) 

Note that in spherical nuclei the tensor part arises due to unpaired nucleons only. The scalar 
part does not require nuclear spin, therefore several nucleons can contribute to it; however, 
if all states in a fine-structure doublet are occupied, then the contributions cancel exactly. 

It is interesting to consider the P, T-odd nuclear polarizability of a deformed nucleus. In 
this case we would expect contributions to the EDM from many nuclei (collective polarizabil- 
ity) because in deformed nuclei about A 2 / 3 nucleons are in open shells. This enhancement 
is similar to that of the P, T-odd nuclear quadrupole moment discussed in [Q. 



II. ESTIMATE FOR INDUCED ATOMIC EDM 

Consider an atom with one external electron above closed shells. We are interested in 
the atomic EDM induced by the interaction of atomic electrons with the P, T-odd perturbed 
atomic nucleus polarized by the fields of the external electron (we use the nuclear model 
defined in the previous section). The induced atomic EDM arises in the third order of 
perturbation theory, 

^ (0, 0'\d z \n',n)(n,n'\H E \m / ,m)(m,m / \H B \0 / ,0) 
datom= 2^ (p w p -— + permutations, (16) 

n,n',m,,m' 1^0 + ^0' ~ &n ~ &n>){&0 + -&0' ~ An ~ &m') 

where unprimed states (m, n, 0) denote electronic states, primed states (m', n', 0') denote 
(P, T-odd perturbed) nuclear states, and the electric and magnetic interactions between the 
external electron and the unpaired nucleon are 

H E = ~^d„ R H B = -jpVn ■ (R x a), (17) 

d n and \x n are the electric and magnetic dipole moments of the nucleus, and a is the electron 
Dirac matrix. We use the Dirac interaction because the external electron close to the nucleus 
is relativistic. 
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To measure an EDM, an external electric field is applied to the system which couples 
to the electron electric dipole moment operator d z . Terms in which this operator admixes 
excited states, for example 

y-, (0,0'\HE\n' ,n)(n,n'\d z \n' ,m)(m,n'\HB\0' ,0) 

m,n,n' (^0 + Eqi — E n — E n ')(E + Eqi — E m — E n i) 

have very large (nuclear) energy denominators as here the electric field interacts with the 
excited states of the atom. We can therefore neglect this contribution to the atomic EDM. 
The atomic EDM induced by the interactions (Eq. [IT]) can therefore be presented as 

, (OjW < n ^ (n,0>\H E \m>)(m'\H B \0',0) 1 

n — &n) m , — An' J 

where we have set n! = 0', because d z does not act on the nuclear wavefunctions. Typical 
electron energies are much smaller than nuclear energies, (Eo—E m ) << (Eo>—E m >), therefore 
we have used closure for the electron states, \ m ){ m \ = 1- Substituting Eq. [17| into Eq. 
15 , and using the definition for the P, T-odd nuclear polarizability (Eq. [I]), we obtain 



, ^ 9 e 2 R ^ (0\d z \n)(n\R t R iap \0) (0\d z \n) (n\H PT \0) 



where we have defined an effective Hamiltonian Hpt which describes the interaction of the 
P, T-odd nuclear polarizability with atomic electrons. 

In this approximation, the contribution of the scalar polarizability to the atomic EDM 
is zero (R • [R x a] =0). However, there is no theorem which forbids this contribution. 
Indeed, consider non-relativistic expressions for magnetic and electric fields produced by a 
charged particle with magnetic moment fi = /is, 

nxv 3(n-u)n— u 8tt „. , „ n , . 

where v is the velocity of the particle, n = H/R. Substituting these expressions into the 
interaction energy 

s ■ n ! Att 

P s 6 tk H t E k = 2efif3 s —i— + —5(R)] (22) 

we see that the first, longer-range orbital field in the expression for H indeed does not give 
any contribution to the interaction. However, the spin magnetic field, which decays more 
rapidly with distance, does contribute. The angular structure of the effective interaction 
(Eq. ^) is similar to that between an electron EDM and an atomic electric field (see, 
e.g., 0). However, the interaction (Eq. ^2|) decays with distance very rapidly, and is very 
singular at small distances. This strong singularity is an indication that one should use the 
relativistic form for the magnetic interaction (Eq. |I7|) in order to describe the small distance 
contribution correctly. We showed above that in this case our calculation gives zero for the 
(3 S contribution to the EDM. 
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The only contribution of the nuclear P, T-odd polarizability to the atomic EDM therefore 
arises due to the tensor term. Inserting the expression for (5^ (Eq. |15|) into the above 
expression (Eq. ^) we obtain for the effective Hamiltonian Hpx 



e 2 



Hpt = -7T^ e N^N(g s - 9i)( — e kl pRiRia p [IiI k + hU - -S ik I(I + 1)]. (23) 
Zri CUfs l 6 

In earlier works the value of atomic EDMs of Cs fTjfJ] , Dy |||5|] and Ra induced by the 

(P, T-odd) nuclear magnetic quadrupole moment (MQM) have been calculated. It is there- 
fore useful to compare Hpt with the Hamiltonian Hmqm which describes the interaction of 
a nuclear MQM with atomic electrons, 

Hmqm = ^5 2/(2/ — fj 6 ^^^^^^^ ^ k ^ i ~ + 1)}, (24) 

in order to obtain an estimate for the value of the atomic EDM induced by the nuclear 
polarizability. Here M is the magnetic quadrupole moment which, for / = I + 1/2, is 
M = -&N(g s - 2gz)(2i~ - 1). The ratio of H PT to H M qm is 

H PT _ 1 (g s - gi)ee N ^ 



H M qm R(g s - Ig^J' 

The largest contribution to the atomic EDM arises due to mixing between s and P3/2 
electron states. The radial integral for the magnetic quadrupole converges at R ~ ap/Z 
||, where ap is the Bohr radius. The Hpt matrix element is more singular at the nucleus 
than that of Hmqm- Considering s and P3/2 mixing, the EDM induced by the nuclear 
polarizability contains an extra factor 

S„, s - 7l t 72 ~^ (^f-) 3 - 7 ""> (26) 
71+72-0 ZR N 

where Pat is the nuclear radius, 71 = y/1 — Z 2 a 2 and 7 2 = \/4 — Z 2 a 2 . For Z 2 a 2 << 1, 

3 Z 2 a 2 ap 

has a very weak dependence on the nuclear radius Rn which we use as a cut-off parameter 
in the integration over P. Therefore, the ratio of Hpp to Hmqm is 

_Hpt_ ^ Ze 2 S re y (2g) 
H M qm a B^is 

For heavy atoms this ratio is of the order of 1%. 

We are grateful to O.P. Sushkov and M.Yu. Kuchiev for useful discussions. This work 
was supported by the Australian Research Council. 
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